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We show that a stacking approach to galaxy clusters can improve current limits on decaying dark 
matter by a factor > 5 — 100, with respect to a single source analysis, for all-sky instruments such 
as Fermi-LAT. Based on the largest sample of X-ray-selected galaxy clusters available to date (the 
MCXC meta-catalogue) , we provide all the astrophysical information, in particular the astrophysical 
term for decaying dark matter, required to perform an analysis with current instruments. 
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I. INTRODUCTION 

The nature of dark matter (DM) is one of the major 
questions in both modern astrophysics and fundamental 
physics. In the most popular extensions of the standard 
model, an explicit symmetry provides stable annihilating 
DM candidates. Nevertheless, decaying DM is an equally 
well-motivated scenario; it arises in several consistent ex- 
tensions beyond the standard model. Among those can- 
didates and frameworks, a non-exhaustive list includes: 
gravitinos in R-parity breaking models [l J2| , sterile neu- 
trinos 3, hidden sector gauge bosons J], hidden sec- 
tor particles [5|, right-handed neutrinos/sneutrinos [6|, 
bound states of strongly interacting particles Q ■ To be 
viable in the DM context, these models should induce a 
DM lifetime longer than the age of the universe. 

For the case considered here (see below), the total flux 
expected from decay, in a given direction {I, b) (Galactic 
coordinates) and integrated over the solid angle Ail, is 
given by the product of a particle physics term with an 
astrophysical term D 



candidates above a few hundred GeV [g]). The mass of 
the candidate can span a large range depending on the 
model. For this reason, a signal is generally searched for 
in X-rays [§, [lO] , or in 7-rays [ll| . In the 7-ray regime, 
the most stringent limit on the DM lifetime (that arises 
in the spectral term) comes from the non-detection by 
the Fermi-LAT collaboration |l2| of emission from the 
directions of galaxy clusters [1J| . 

Here we focus on the astrophysical factor 
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Depending on the decaying DM candidate, the photon 
spectrum dN / dE can be a mono-energetic line (from 77 
and ^v channels), or a continuum (e.g., for some gravitino 
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which is the integral of the DM density, p{l' , fi), over line 
of sight V and solid angle Afi = 27r- (1 — cos(aint)), where 
aint is the integration angle. When D is computed over 
cosmological distances, the spatial term becomes cou- 
pled to the energy-dependent term as 7-rays are absorbed 
alon g th e line of sight (e~'^(-^'^) attenuation factor, see, 
c-g-, uM)- 'Cosmological' 7-rays will also be redshifted, 
affecting the spectrum. Here however, our analysis relies 
on the MCXC catalogue of galaxy clusters [1^ , the red- 
shift distribution of which peaks at 2 ~ 0.1 (see their Fig. 
1) so that we can safely neglect the above processes. This 
allows us to factor out energy-dependent effects, i.e. the 
spectrum resulting from a prompt emission and also an 
inverse Compton (IC) contribution from scattered DM- 
induced electrons and positrons. The benefit is twofolds: 
i) it simplifies the discussion, and ii) also avoids intro- 
ducing strongly DM model-dependent factors. Note how- 
ever, that the respective IC contributions in the Galaxy 
and in the clusters are probably different, and this should 
be kept in mind when comparing the Galactic DM diffuse 



and galaxy cluster D-iactors. 

As in the case of searches for annihilating DM, in- 
teresting targets include dwarf spheroidal galaxies, our 
own Galaxy, and clusters of galaxies [ia| • Recently, dif- 
fuse 7-ray emission |14l | and cross-constraints with other 
channels (e.g^ anti-protons and positrons) have also been 
considered |17H21| . These studies rely on the analysis of 
single objects or small samples constituted by what are 
thought to be the best targets. As underlined in Q, one 
advantage of combining several galaxy clusters is that 
uncertainties (for example on the total mass and DM 
concentration) somewhat cancel out. 

With this in mind, we show below how the astrophys- 
ical term D can be increased by a stacking strategy of 
galaxy clusters without compromising sensitivity due to 
increasing the amount of background integrated. For 
that purpose, we take advantage of the recently pub- 
lished Meta-Catalog of X-ray detected Clusters, MCXC 
[15l |. which encompasses 1743 clusters of galaxies (mak- 
ing it the largest sample of clusters with detected X-ray 
emission from their hot gas). 
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FIG. 1: Astrophysical factor for decaying DM. Symbols are 
for galaxy clusters, solid line is for the Galactic DM back- 
ground, and the dotted and dashed lines are an upper-limit 
(no absorption on the extragalactic background light) for the 
isotropic extragalactic DM background (see Eq. [3]). The 10 
brightest objects are identified with coloured symbols. 



10' 



lo- 



ll. METHOD 

A. Galactic and cluster signal 

All calculations of the Z?-factor are performed with the 
public code CLUMPY v20n.09 Q. The Galactic DM 
halo is chosen to follow an Einasto profile [2J| with a local 
density p© = 0.3 GeV cm~^. We assume that the DM 
profiles of all the galaxy clusters in the catalogue follow 
a NFW paramctrisation |2J]. For a cluster, we use the 
standard definition for Ra to be the radius within which 
the average density reaches A times the critical density 
of the Universe (at a given redshift). The mass Ma is 
then simply the mass enclosed within Ra- Most observa- 
tional constraints and predictions are expressed in terms 
of A = 500 or A = 200. Using values for M500 (de- 
rived from the X-ray luminosity) provided in the MCXC 
meta-catalogue |15| . we use a mass-concentration rela- 
tionship [25| to derive the scale radius and normalisation 
of the NFW profile for each cluster. To first order, the 
uncertainties related to X-ray observations (i.e. ~ 15% 
average on M500), and the scatter related to the mass- 
concentration relationship [25| average out if the stacked 
signal from a large sample of clusters is considered. 



B. Extragalactic isotropic signal 

As shown in [ij], contrary to the DM annihilation 
case, the extra-galactic DM decay contribution is a ro- 
bust quantity depending mostly on the mean DM density 
of the universe r^DMPc.o- For the sake of simplicity, we 
neglect the absorption t(E,z) as it has only a moderate 
impact below 100 GcV [ij]. Doing so yields an upper 



limit for this contribution that is written as 
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where H is the Hubble constant for the concordant A- 
CDM cosmology. This quantity is used below for illus- 
tration purposes only, as a limit to the continuum signal 
which can drown the cluster signal (see in particular Sec- 
tion HUB]). 



III. RESULTS 

Figure [1] shows the astrophysical term D (in 
Mq kpc~^) for each cluster (plus symbols) as a function 
of its angle 9 w.r.t. the Galactic centre direction. A large 
number of objects outshine the Galactic DM background 
(the solid line), fewer outshine the extragalactic signal 
upper limit (dotted line). Being a compilation of various 
X-ray cluster catalogues, the MCXC meta-catalogue is 
not complete at any redshift or in any mass range. This 
diversity is reflected by the less populated region below 
the Galactic signal (solid line) in Fig. [T] In principle, 
an integration along z > of the extragalactic signal 
amounts to some double counting of the MCXC galaxy 
clusters. The dashed line calculated for z > 1 illustrates 
how much this double counting could be at most (since 
the MCXC catalogue is not complete up to this redshift). 

An integration angle of aint = 0.1° is adopted, cor- 
responding to the typical angular resolution of current 
7-ray instruments well above threshold (e.g., Fermi-LAT 
for energies above - 10 GeV, H.E.S.S. above -300 GeV). 
This choice is appropriate for the signal from cluster ha- 
los, as their typical angular scale (i.e., i?5oo) is 0.1° — 1°. 
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FIG. 2; logA'^ — logD for the galaxy cluster population for 
aint = 0.1°. The behaviour for the brightest objects is con- 



sistent with N (X D 
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However, we underline that the Galaetic and extragalac- 
tic signals scale with a?^^ , so that smaller integration an- 
gles are in principle favoured to increase the contrast be- 
tween the cluster signal and these DM backgrounds. 



A. log iV - log D 

The promise of a stacking analysis for decay is sug- 
gested from Fig. [2l The log N — log D plot has a slope 
of —2.3 at large values of D: the number of objects in- 
creases faster than the signal decreases. The incomplete- 
ness of the MCXC meta-catalogue is seen as a drop for 
D < 10^ Mq kpc~^. A larger and more complete cata- 
logue with a well defined selection function could further 
increase the gain provided by a stacking approach if the 
A^ ex D^^-^ behaviour continues down to smaller D val- 
ues. In a few years from now, the eROSITA mission [2g| 
should provide such a catalogue. There is naturally a 
limit to this gain, as at large redshift, clusters may not 
have formed yet. The maximum gain can in principle 
be estimated from the theoretical redshift distribution 
of clusters, or from the use of more complete, e.g. op- 
tical, galaxy cluster catalogues (although these observa- 
tions poorly constrain the cluster DM content). For the 
current study wc focus on X-ray identified clusters to re- 
duce uncertainties in the DM signal and to provide a list 
of target objects to current observatories. 



B. Stacking 

Figure [3] (upper panel) shows the cumulated signal of 
all galaxy clusters having an astrophysical factor > D a,s 
a function of D. In all cases the cumulative signal satu- 
rates at < 800 clusters, consistent with the drop-out seen 
in Fig. [5J For the case of aint = 0.5° (black circles) the 
integrated D for 791 clusters is 6-10® Mq kpc~^, a factor 
^ 30 above the value for the top two clusters (furthest 
right point). For each D bin, the corresponding cumula- 
tive background signal from the Galactic DM halo (resp. 
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FIG. 3: Upper panel: cumulative galaxy cluster signal (solid 
symbols) and Galactic (solid lines) and extra-galactic upper 
limit (dashed lines) DM background as a function of D. The 
three sets of curves correspond to three different integration 
angles: aint = 0.5° (black circles), aint = 0.1° (blue squares), 
and Qint = 0.05° (red stars). The number above a given 
bin corresponds to the total number of clusters with a signal 
larger than D. The far left-hand bins contain a total of 1743 
clusters. Bottom panel: cumulative D divided by ainfs/iV, 
proportional to the signal-to-noise ratio for a given cluster 
sample. 



extragalactic) is also shown with solid lines (resp. dashed 
lines). The region between the two lines provides the 
range of accumulated diffuse DM signal to expect. For 
29 clusters, the accumulated background signal from the 
Galactic DM halo reaches the level of the cluster signal. 
The integration angle of 0.5°, as used in Fermi-LAT anal- 
yses of dwarf spheroidal Galaxies [23, ^M j corresponds to 
the Fermi-LAT angular resolution at around 1 GeV, de- 
grading to ~ 1° at lower energies and reaching 0.1° above 
10 GeV. The blue squares and red stars show the result 
of repeating the analysis for ai^t ~ 0.1° and 0.05° re- 
spectively (a resolution generally within the capabilities 
of X-ray instruments, or within reach of the next genera- 
tion of air-Cerenkov 7-ray telescopes). In these cases the 



accumulated Galactic DM plus extra-galactic signal is al- 
most always below the cumulative cluster signal. Because 
of a smaller integration angle, saturation is not reached 
at the same cumulated value. For aint = 0.1° a gain of a 
factor ^ 100 in signal is observed between stacking two 
and 948 sources in the cumulative D. 

The best observational strategy depends not only 
on the integrated signal and Galactic halo plus extra- 
galactic DM background, but on the (generally much 
higher) level of charged particle and diffuse astrophys- 
ical gamma-ray backgrounds. In the background-limited 
regime the best approach is to maximise the signal- 
to-noise ratio or signal divided by the square-root of 
background. In the case of uniform background the 
signal-to-noise is proportional to the cumulative D di- 
vided by aintV^V (where N is the number of clusters 
with Di > D), as shown in the bottom panel of fig. [3l 
The peaks in these curves (at 29, 948 and 971 clusters 
for aint = 0.5°, 0.1° and 0.05° respectively) correspond 
to the optimum stacking approach in the background- 
limited, rather than signal-limited, regime. With the cur- 
rent MCXC meta-catalogue, there seems to be no advan- 
tage in resolutions better than 0.1°. In the signal-limited 
regime a larger integration angle is preferred. Stacking 
leads to a signal-to- noise increase of a factor ^ 5. With a 
more complete galaxy cluster catalogue, the 'saturation' 
regime (the break seen in the upper-panel curves) would 
be moved towards smaller D values and the benefit of 
the stacking approach could be greater. 



IV. DISCUSSION AND CONCLUSION 

As well as (close to uniform) charged particle back- 
grounds and the diffuse Galactic DM and astrophysical 
backgrounds (see e.g. [2al)j non-DM gamma-ray emission 
from galaxy clusters must be considered as another pos- 
sible source of background and/or confusion [see 29, '3^ . 
No individual galaxy cluster ~ but for Virgo, for which a 
gamma-ray emission was detected positionally consistent 
with the centre of M87 by Fermi-LAT [U - has been 
detected in gamma-rays so far, but accelerated particles 
are expected to provide an additional gamma-ray signal 
in clusters. Given the present situation of non-detection 
in the direction of individual massive halos, a stacking 
analysis is advantageous to place limits on DM decay. 
Any object for which a signal (of suspected non-DM ori- 
gin) is observed can be discarded in the analysis at little 
cost for the exclusion limit. On the other hand, if any 
signal is observed and is consistent with an astrophysical 
origin, stacking or merely looking at galaxy clusters will 
be extremely challenging for DM searches. Alternatively, 
the strong scaling log N — log D could be a good diagnos- 
tic to disentangle a decay from an astrophysical signal if 
the latter shows a different scaling (Maurin et al., sub- 
mitted). The case of annihilating DM, which is found to 
be less favourable but more crucial physically (limits of 
non-detection are close to the cross-section expected from 



that required to match the thermal relic abundance) is 
discussed in Nezri et al. (in preparation), focusing on the 
instrumental response of current and future observatories 
(Fermi and CTA). 

To ease further analysis based on our model, we pro- 
vide with this paper a file containing all the necessary 
inputs for each cluster. An analysis on real data can 
certainly be optimised by adapting the integration re- 
gion for each cluster. Assuming all galaxy clusters share 
the same DM profile, and given the mass range spanned 
by the MCXC, we can make the first order approxima- 
tion that their concentration parameters are the same. 
In other words, the concentration for an NFW profile 
is c{M) = Rvir/rs and we assume c(lO^*M0) ~ 5 [2^. 
Defining 
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we find that there is a universal dependence of the frac- 
tion 
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on the DM-decay signal. Hence, given £'(0.1)° and as 
(available in the supplementary file |37| . a sample of 
which being given in table U) , and a parametrisation of 
this universal dependence 
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the term D can be calculated for any integration angle, 
using 
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This parametrisation describing the fraction of the signal 
in a given angular region is valid down to J-d = 10~^. 
We note that this function is unchanged if an Einasto 
profile is used, but the D- values obtained from such a 
profile are larger than that calculated with an NFW DM 
distribution. 

The values for the ten brightest clusters are given in 
Table HI with values of D provided for a 0.1° integra- 
tion angle (see Supplemental Material [37] for the full 
list of clusters). We have repeated (not shown) the cal- 
culation for aint = 1° (an angle for which most clusters 
fall below the Galactic DM background) to compare the 
D factors with those of Q, UM'- o^r values for Fornax, 
Coma, AWM 7 and NGC 4636 are found to be lower, but 
within a factor 2 of what was found in these previous 
works. In particular, the large difference for the Fornax 
cluster explains why these authors have flagged it as be- 
ing the best target for DM decay when it does not make 



TABLE I: Ten brightest galaxy clusters (positions taken from 
the MCXC meta-catalogue ^) in DM-decay for aint = 0.1°. 



Name 


Index 


I 


b 


d 


Os 


D(0.1°) 




(MCXC) 


(deg) (deg) (Mpc) 


(deg) (Mq kpc-2) 


A426 


258 


150.6 


-13.3 


75.0 


0.44 


1.0- 10* 


Virgo 


884 


283.8 


74.4 


15.4 


1.09 


1.0 


10* 


Coma 


943 


57.2 


88.0 


96.2 


0.29 


7.3 


10^ 


Ophiuchus 1304 


0.6 


9.3 


116. 


0.27 


7.3 


10^ 


A3526 


915 


302.4 


21.6 


48.1 


0.39 


6.8 


10^ 


A3627 


1231 


325.3 


-7.1 


66.0 


0.32 


6.5 


10^ 


AWM7 


224 


146.3 


-15.6 


72.1 


0.28 


5.7 


10^ 


A1367 


792 


235.1 


73.0 


89.3 


0.24 


5.4 


10^ 


A3571 


1048 


316.3 


28.6 


160. 


0.18 


5.4 


10^ 


A2199 


1249 


62.9 


43.7 


124. 


0.20 


5.1 


10^ 



it in our top-ten. This is understood as follows. De- 
cay is very sensitive to the mass estimate (which goes 
directly as the integration of the density) . In |^, [l^ ^ the 
authors use Afsoo values from the HIFLUGCS catalogue 
[33 . [33j that are larger than the ones provided in the 
MCXC catalogue (e.g., factor of ^5, 2, and 2 for For- 
nax, Coma, and AMW7 respectively). For the former 
analyseSjthe authors used a /3-model for the gas distri- 
bution [32|, |33], whereas the MCXC relies on the more 
accurate AB-model [SJI. As discussed in the App. A of 
[15| . beta-models can produce factor 2 higher, or lower, 
values of M500, depending on the core radii used. The 
masses for the MCXC sample could still suffer from sys- 
tematics, as for instance from the use of the hydrostatic 
equilibrium hypothesis. However many comparisons to 
numerical simulations indicate that these uncertainties 



are < 15 — 20% [3a]. In any case, if the signal from the 
brightest source is overestimated (resp. underestimated) 
by a factor / the improvement factor will increase (resp. 
decrease) roughly by the same amount. For the MCXC, 
we recall that M500 and -R500 are obtained from the X-ray 
luminosity, making use of the Lx — Mnnn scaling relation 
derived for the REXCESS sample by [3^|. The mass esti- 
mate not only links to the Lx measurement, but also to 
the intrinsic scatter about the scaling relation. Further 
conversions to i?200 niay amplify this dispersion effect. 
This stresses that caution is advised when working with 
individual clusters for which mass and radius are inferred 
from the X-ray luminosity. In that respect, stacking pro- 
vides a more robust approach as it washes out those un- 
certainties. 

To conclude, we have shown that a stacking ap- 
proach can help to push down by a factor of 5 (for 
the background-limited regime) to 100 (in the signal- 
limited case) the current observational limits on decaying 
DM (when using the MCXC meta-catalogue of galaxy 
clusters). A more thorough analysis (e.g., using deeper 
galaxy cluster catalogues and including absorption on the 
extragalactic background, uncertainties on DM profiles, 
etc.) is left for future work. Meanwhile, we provide the 
necessary inputs to apply this idea to existing experi- 
ments (for example Fermi-LAT). 
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